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Experimental 


A  considerable  amount  of  nev  data,  including  x-wire  measurements, 
has  been  obtained  in  the  turbulent  vake  behind  a  flat  plate.  Sinous 
disturbances  at  several  amplitudes  and  frequencies  were  introduced  into  the 
vake  by  oscillating  a  small  trailing  edge  flap.  The  Strouhal  numbers  of 
the  perturbations  were  specially  chosen  so  that  the  down  stream  location  of 
the  neutral  point  (where  the  spatial  amplification  rate  obtained  from 
linear  stability  theory  vanishes)  vas  veil  within  the  the  range  of 
measurements.  The  streamvise  variation  of  the  half  width  of  the  vake,  L0» 
and  the  centerline  velocity  deficit,  u0,  were  dependent  on  the  amplitude 
level  for  a  fixed  frequency.  Large  deviations  from  the  well  known  square 
root  behavior  are  observed  at  large  amplitude  forcing  levels  starting  in 
the  region  of  the  neutral  point.  This  is  shown  in  Figure  1  where  L0 
<Wuo0  ,a  nondimen sional  growth  variable  is  plotted  versus  (x-xo)/20  for  a 
forcing  frequency,  f,  of  70  hr  and  several  forcing  amplitudes.  U^  is  the 
free  stream  velocity,  ©  is  the  momentum  thickness,  and  xQ  £s  the  virtual 
origin.  ©  remains  effectively  constant  for  these  cases.  The  unforced 
case  shows  linear  behavior  as  the  ratio  of  L0/U{)  £8  proportional  to  x,  the 
downstream  coordinate.  The  rate  of  growth  for  the  forced  cases  is 
initially  larger  than  that  for  the  unforced  case.  Depending  on  the  level 
of  forcing,  however,  the  vake  stops  growing  at  some  downstream  location  and 


8?  4 


appears  to  be  in  a  state  of  nearly  parallel  flow.  In  the  high  amplitude 
case  the  wake  even  contracted  some  and  then  resumed  growing  at  a  slower 
rate.  The  downstream  location  where  the  wake  stops  growing  in  each  case 
roughly  corresponds  to  the  location  of  the  neutral  point  for  that  case.  The 
exact  location  of  the  neutral  point  is  difficult  to  determine  as  non-linear 
terms  are  not  negligible  near  the  end  of  the  amplified  region.  The 
variation  in  the  growth  rate  of  the  wake  near  the  neutral  point  is  caused 
by  the  nonlinear  interaction  of  the  perturbation  wave  with  the  mean  flow. 
In  the  region  upstream  of  the  neutral  point,  the  disturbances  are  being 
amplified  and  energy  is  transferred  from  the  mean  flow  to  the  disturbance. 
The  wave  induced  stress  changes  sign  through  the  neutral  point,  so 
downstream  of  the  neutral  point  energy  is  transferred  from  the  disturbance 
to  the  mean  flow.  This  demonstrates  that  the  lateral  rate  of  spread  of  the 
wake  is  closely  linked  with  the  growth  of  the  disturbance  wave; 
amplification  of  the  wave  results  in  a  transfer  of  energy  from  the  mean 


flow  and  a  divergence  of  the  mean  flow.  Further  evidence  to  substantiate 

this  can  be  obtained  from  noting  that  the  initial  divergence  of  the  wake 

increases  with  increasing  forcing  amplitude.  The  contraction  of  the  wake 

may  be  attributed  to  the  transfer  of  energy  from  the  disturbance  wave  to 

the  mean  flow  in  the  damped  region  beyond  the  neutral  point. 

To  consider  some  of  the  interesting  features  of  the  downstream 

development  of  a  sinuous  wave,  let's  examine  the  case  of  a  medium  amplitude 

80  hz  wave.  Figure  2  shows  the  nondimen sional  growth  variable  /u0  0 

0 

versus  (x-xQ)/2  .  The  measured  mean  velocity  profiles  for  various 
downstream  locations  are  presented  in  Figure  3  in  self  preserving  form. 
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Figure  3a  shows  the  profiles  for  the  downstream  region  labelled  a.  These 
profiles  agree  well  with  those  for  the  unforced  case  shown  by  the  solid 
curve,  which  represents  a  curve  fit  to  the  unforced  profiles.  The  profiles 
for  the  region  b,  located  just  past  the  neutral  point,  differ  from  the 
unforced  data  as  shown  in  Figure  3b.  Further  downstream  in  region  c,  the 
changes  are  more  significant.  The  three  regions  defined  in  Figure  3 
correspond  to  the  amplified,  neutral  and  damped  regions  respectively  in  the 
context  of  linear  stability  theory. 

Solution  to  the  inviscid  Orr-Summerf eld  equation  using  the  appropriate 
measured  mean  velocity  profile  gives  the  eignefunction  distributions 
typical  for  each  region.  In  Figure  4  these  solutions  are  presented  in  the 
form  of  the  amplitude  and  phase  of  the  u  component  of  the  disturbance  wave. 
Note  the  dramatic  change  in  the  phase  of  the  u  component  in  the  damped 
region..  The  phase  of  the  v  component  does  not  change  sign,  and  therefore 
the  sign  of  the  wave  induced  or  coherent  Reynolds  stress  changes  in  the 
neutral  region.  The  computed  wave  induced  or  coherent  stress  distributions 
are  shown  in  Figure  5.  Four  plots  are  presented  for  the  neutral  region 
showing  the  downstream  evolution  of  the  change  in  sign.  In  Figure  6, 
typical  measured  amplitude  and  phase  distributions  of  the  u  component  for 
each  region  are  presented  for  comparison.  The  measured  coherent,  turbulent 
and  total  Reynolds  stress  for  the  three  regions  are  displayed  in  Figure  7. 
The  vertical  markers  above  the  data  presented  in  Figure  2  indicate  the 
measurement  locations  for  the  data  shown.  Three  data  sets  are  shown  for 
the  neutral  region  demonstrating  the  change  in  sign  of  the  coherent  stress. 
Each  distribution  has  been  normalised  by  the  maximum  value  of  the  total 
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stress  corresponding  to  that  downstream  location,  so  the  relative  amount  of 
coherent  and  turbulent  stress  can  be  determined.  Figure  7a  shovs  the 
measured  stresses  in  the  amplified  region.  Notice  the  sign  of  the  coherent 
stress  is  the  same  as  that  of  the  turbulent  stress,  while  its  magnitude  is 
roughly  twice  that  of  the  turbulent  stress.  In  Figure  7e,  the  measured 
stresses  in  the  damped  region  c  are  presented.  Note  that  sign  of  the 
coherent  stress  has  changed  and  its  magnitude  is  now  only  one-third  of  that 
of  the  turbulent  stress. 

Figure  8  presents  the  downstream  development  of  the  integral 
amplitudes  of  u^,  v^,  and  uv  in  the  wave  component  at  the  fundamental 
frequency.  Each  variable  is  decomposed  into  its  coherent  (phase  averaged) 
and  turbulent  parts  and  non-dimensional ized  by  uq  and  L0*  ^h«  peak  in  the 
v  component  of  the  fundamental  occurs  at  about  the  downstream  location  of 
the  neutral  point  and  corresponds  to  the  position  where  the  growth  rate  of 
the  wake  changes  as  can  be  seen  from  Figure  2.  Both  u  and  uv  peak  sooner 
than  v,  for  both  the  fundamental  and  harmonic.  Note  that  at  large  x,  the 
coherent  stress  is  vanishing  while  the  turbulent  part  is  increasing  and 
approaching  a  value  which  is  characteristic  of  the  self  preserving  unforced 
wake.  A  comparison  of  the  evolution  of  the  fundamental  and  harmonic  for  v 
is  shown  in  Figure  9.  As  is  readily  apparent,  the  harmonic  is  more  than  an 
order  of  magnitude  smaller.  The  measured  amplitude  and  phase  distributions 
of  the  U  component  of  the  harmonic  are  shown  in  Figure  10. 


Numerical  Simulations 


In  the  past  year,  ve  made  considerable  progress  towards  development  of 
a  numerical  method  that  will  be  applicable  for  investigations  of  the 
dynamics  of  coherent  structures  in  turbulent  vakes.  As  mentioned  in  last 
year's  report  for  the  actual  development  of  the  numerical  methods  for 
solving  the  complete  Navier-Stokes  equations,  we  decided  to  use  a 
different,  but  related,  flow  geometry,  namely  that  of  a  free  shear  layer. 
The  main  motivation  for  this  was  the  fact  that  for  free  shear  layers  a 
considerably  larger  body  of  substantiated  experimental  results  are 
available  than  for  wake  flows.  Therefore,  a  thorough  check  of  our  new 
numerical  methods  is  possible  before  applying  it  to  investigations  of  wake 
flows.  The  development  of  a  numerical  method  which  was  based  on  a 
"temporal"  approach  (disturbances  simplify  or  decay  in  time)  was  already 
concluded  some  time  go,  and  we  reported  it  in  last  year's  report.  During 
the  past  year,  we  concluded  the  development  of  a  numerica*l  method  for 
solving  the  complete  Navier-Stokes  equations  which  is  based  on  the 
"spatial"  approach.  (Fh.D.  Thesis  by  C.  D.  Pruett)  The  method  enables 
realistic  simulations  of  linear  and  nonlinear  development  of  stability 
waves  where  development  (growth  or  decay)  of  the  disturbances  is  spatial 
(in  downstream  and/or  spanwise  direction)  as  in  real  laboratory 
experimental  situations.  In  particular,  since  the  nonlinear  effects  are 
fully  incorporated,  it  allows  us  to  investigate  the  saturation  of  the 
instability  waves,  and  therefore  the  intensity  of  the  coherent  structures 
can  be  predicted.  In  addition,  due  to  the  three-dimensional  capability,  we 
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can  also  investigate  the  importance  of  three-dimensional  effects  or  study 
the  breakdown  of  two-dimensional  structures  into  three-dimensional  ones. 

Typical  results  obtained  from  simulations  of  a  transitional  shear 
layer  as  shown  in  Figure  11  and  Figure  12).  Figure  11  shows  the  downstream 
growth  of  various  wave  components  for  a  forced  transitional  free  shear 
layer.  Comparison  with  experiments  by  Hiksad  (1972)  shows  excellent 
qualitative  and  quantitative  agreement.  Also,  for  the  transitional  free 
shear  layer  Figure  12  exhibits  the  spatial  development  of  the  three- 
dimensionality  (development  of  secondary  instability). 

We  carried  out  other  comparison  calculations  with  other  shear  layer 
experiments  to  convince  ourselves  that  our  numerical  method  works  properly. 
So  far,  we  have  had  only  positive  experiences  and,  therefore,  we  conclude 
that  we  are  prepared  to  apply  our  method  now  to  wake  simulations.  In 
summary,  together  with  the  numerical  methods  we  have  developed  and  tried 
previously,  we  have  the  following  computer  codes  that  can  be  readily 
applied  to  wake  investigations: 

a.  Two-Dimensional,  Temporal  (Spatially  Periodic) 

Finite  Difference  in  y.  Finite  Difference  in  x 
Finite  Differences  in  y,  Pseudo-Spectral  in  x 

b.  Two-Dimensional,  Spatial  (Amplification  in  x) 

Finite  Differences  in  x  and  y 

c.  Three-Dimensional,  Temporal  (Periodic  in  x) 

Finite  Differences  in  y  and  x,  Pseudo-Spectral  in  z  (Spanwise) 
Finite  Differences  in  y,  PrTidj-S* jctral  in  x  and  z 


F  =  7 0 H Z  VARIOUS  FORCING  LEVELS 


Figure  1.  The  downstream  development  of  the  growth  parameter  for 
and  various  forcing  levels. 


Figure  2.  The  downstream  evolution  of  the  growth  parameter  for  a  medium 
amplitude  forcing  and  f  a  00  itz.  Downstream  regions  defined. 


Fi(|iire  4.  Computed  amplitude  and  phase  distributions  for  u  component  of  disturbance  wave 
typical  of  amplified,  neutral  and  damped  regions  from  linear  stability  theory. 
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figure  7.  Typical  measured  distributions  of  total,  wave-induced  or  coherent,  and  turbulent  Reynolds 
stresses  for  the  downstream  location'',  identified  in  Figure  2. 
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Figure  11.  Growth  of  U|-ve1ocity  spectra  maxima  (Inset  from  Mlksad,  1972) 
Transitional  free  shear  layer 


+2 


Contour  interval  «  0.0015  (zero  controu  suppressed) 


Figure  12.  Evolution  of  strearrwise  vorticity  a)  n»50,  b)  n*58, 
c)  n»66,  d)  n»86. 

(n»x)/ix)  Transitional  free  shear  layer. 


